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(54) P«>c6d6 de f lltrage spatial du bruit dans une Image num6rlque, et disposltif mettant en 
oeuvre ce proc6d6 

(57) L'invention concerne un proc6d6 de filtrage 
spatial du bruit dans une image num6rique. Ce proc6d6 
comprend la localisation d'un pixel couranl tPo(xo. yo)l 
d-intensit6 brurt6e [lo(xo.yo)l : 'a d6finition d'une fen§tre 
glissante [WD(xo. yo)l centr6e sur le pixel courant ; la 
selection d'un nombre (N). fixe ou al6atoire, de f.ltres 
6l6mentaires al6atoires orient6s [Fi. Fa.... Fn) ayant un 

axe d'orientation IPqZi. P0Z2 PoZnI Passant par ce 

pixel courant, et munis de caract6ristiques dont au 
moins une est alfeatoire. incluant : un angle d'orientation 

[0^. ©2 ©n1 entre I'axe d'orientation et un axe de 

r6f6rence (OX), et une 6tendue de support d6finie en 
nombre de pixels ; la combinaison des intensity des 
pixels sur les supports des filtres 6l6mentaires indivi- 
duellement pour lournir des intensit6s 6l6mentaires ; et 
la combinaison des intensit6s 6l6mentaires pour fournir 
une intensity filtrfee [Ro(xo. Yo)! ^ appliquer au pixel 
courant. 

L'invention concerne aussi un dispositif de traite- 
ment d'image mettant en oeuvre ce proc6d6. 
Application : imagerie m§dicale num6ris6e 
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Description 

L'invention cx)ncerne un proc6d§ de filtrage spatial 
du bruit dans une image num^rique. 

L'invention concerne 6galement un dtspositif met- 5 
tant en oeuvre ce proc6d6. 

L'invention trouve particuli^rement son application 
dans le trartement d'images m^dicales form§es dans un 
mode fluoroscopique de rayons X. au moyen d un sys- 
teme ou le falsceau de rayons X a une faible intensity. io 
fournissant par suite des images bruit^es et peu con- 
trast6es qui doivent §tre filtr^es afin d en extraire le 
bruit, sans en d^t^riorer les details. 

Le mode fluoroscopique de rayons X est utllis6 par 
exemple pour suivre en temps r§el une operation m^i- is 
cale durant laquelle un outil de diam^tre extremement 
petit tel qu'un cathdter est introduit ou diplac^ dans la 

Un proc§d6 de filtrage du bruit dans une image 
num^rique est 66\k connu par la publication intitul6e "A 20 
new class of Detail-Preserving Filters for Image Proces- 
sing" par Ari NIEMINEN. et alii, in "IEEE TRANSAC- 
TIONS ON PATTERN ANALYSIS AND MACHINE 
INTELLIGENCE VOLPAMl-Q. NM, JANUARY 1987. 
PP 74-90". 25 

Cette publication d^crit un proc6d6 de filtrage utili- 
sant une nouvelle classe de filtres spatials. Cette nou- 
velle classe comprend des structures de filtres Iin6aires 
^ r^ponse impulsionnelle, appel§s "FIR". Le proc6d6 
divutgue par ce document cit6 inclut des combinaisons 30 
'des resultats des filtres lin^aires dans des filtres 
medians, avec plusieurs niveaux de combinaison. II est 
enseigne qu'un tel proc6d6 permet de preserver les 
bords des objets dans I'image. et ainsi de preserver les 
details. Ce proc6d6 inclut la d^inition d'une fenetre glis- 35 
sante sur I'image num^rique k traiter centr6e sur un 
pixel courant, et la definition de masques dans cette 
fenetre pour permettre la realisation d'un certain nom- 
bre de moyennes des intensites des pixels situ6s autour 
du pixel central de la fenetre. Par exemple. une 4o 
moyenne est r^alis^e sur les intensites des pixels situ^s 
sur la m§me llgne de pixels que le pixel central, une 
moyenne est realis^e sur les Intensites des pixels situes 
sur la meme colonne que le pixel central ; les moyennes 
d'intensite peuvent aussi §tre realisees sur des pixels 45 
des diagonales de la fen§tre, cette derniere etant carree 
dans Texemple decrit. Les resultats de ces moyennes. 
qui constituent des filtrages lineaires (FIR) appeies fil- 
trages de base, sont combines dans des filtres medians 
(FMH). Par exemple, les moyennes calcuiees sur les so 
intensites des pixels de la ligne du pixel central k gau- 
che et k droite du pixel central respectivement, sont 
combinees dans un f iltre median avec Tintensite du pixel 
central, et un premier resultat est fourni .; puis les 
moyennes calcuiees sur les pixels de la colonne du 55 
pixel central vers le haut et vers le bas respectivement 
sont combinees avec I'intensite du pixel central dans un 
filtre median, et un second resultat est fourni ; puis les 
premier et second resultats sont combines dans un filtre 



median avec I'intensite du pixel central, et un autre 
resultat est fourni. Seion le nombre de filtrages de base 
calcuies dans la fen§tre, et le nombre d'entrees des fil- 
tres medians, le nombre d'etapes de combinaisons pour 
obtenir le resultat, appeiees niveaux de combinaisons, 
peut §tre plus ou moins grand. 

Un autre proc6de de filtrage du bruit dans une 
image numerique est egalement dej^ connu par la 
publication intituiee "Performance Evaluation of Some 
Noise Reduction Methods" par Wen-Yen WU et alii, 
dans "GRAPHICAL MODELS AND IMAGE PROCES- 
SINGS. vol.54. N°2. March 1992. pp.134-146. 

Cette seconde publication decrit differents proce- 
des de filtrage du bruit incluant d'une maniere generate 
d'abord la selection, selon des criteres donnes, de 
pixels dans une image numerique sur I'intensrte des- 
quelo on r^eliso uno combinaiaon. Qttt« combtnalson 
ptmm mnn uhft ^i^Pti^thnii^h lir^#air0 imWB «|U*uiia 
moyenne. A cet effet. ce procede inclut la definition 
d'une fendtre glissante sur Timage numerique k traiter. 
centree sur un pixel courant, et la selection de pixels 
dans cette fenetre. Un critere de selection peut porter 
sur rintensite des pixels seiedionnes. Dans un exemple 
(SIG), les pixels seiectionnes sont ceux dont la diffe- 
rence d'intensite vis-^-vis du pixel courant situe au cen- 
tre de la fen^re. est inferieure k un seuil donne qui 
depend de la variance du bruit (o^) ; dans un autre 
exemple (KNN), un nombre K donne de pixels ayant les 
intensites les plus proches de cellesdu pixel courant est 
seiectionne. Apres I'operation de selection des pixels 
dans la fenetre glissante carree, leur intensite est com- 
binee soit en effectuant une moyenne (MEAN), soit par 
un filtrage MEDIAN (MED). 

Un probieme qui se pose dans le traitement du bruit 
des images en fluoroscopie est que ces images con- 
tiennent en general des details tres fins ou la represen- 
tation d'objets tres petits, comme par exenple un 
catheter introduit lors de I'operation par le praticlen, ce 
catheter ayant la forme d un fil dont la dimension trans- 
versale est de quelques pixels. Dans ce domaine, II est 
done indispensable d'avoir recours k une methode de 
filtrage du bruit qui conserve les details. 

Toutes les methodes decrites dans les deux publi- 
cations citees prec6demment au titre d'etat de la techni- 
que sont capabtes de fournir un filtrage de Tlmage 
bruitee en conservant les details. 

C^endant, les methodes connues des documents 
cites posent toutes un m§me probieme ; les images 
reconstruites k Taide de ces methodes montrent des 
TACHES (PATCHES). Ceci est dO au fait que ces 
methodes reposent sur I'utilisation d'une fen§tre glis- 
sante dont le centre est deplace d'un pixel sur I'image 
apres la mise en oeuvre de chaque ensemble d'etapes 
de calcul decrit dans Tune ou I'autre des publications 
citees. Ainsi. Taire de recouvrement des fenetres cor- 
respondant k deux pixels adjacents est importante, d'ou 
il resulte que la contribution du bruit residuel dans ces 
deux fenetres est tres correiee. 

Un but de la presente invention est de proposer un 
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proc^6 de filtrage spatial qui ne g6n6re pas de 
TACHES de bruit. 

Un but particulier de rinvention est de proposer un 
tel proc6d6 capable de produire une Elimination tr^s 
pouss6e du bruit, avec preservation des contours des s 
objets dans I'image, et avec preservation du contraste 
derimage. 

Ces buts sont atteints par un proc6d6 de filtrage 
spatial du bruit dans une image ayant la forme d une 
matrice bidimensionnelle de pixels ayant des valeurs io 
d'intensite num6ris6es bruit6es, ce proc6de compre- 
nant des etapes de : 

a) selection, dans une fenfitre glissante centr^e sur 

un pixel courant d'un nombre (N 0) de f iltres 616- 75 
mentaires orient6s ateatoires ayant des supports 
avec un axe d'orientation passant par le pixel cou- 
rant et des caracteristiques dont une ou plusieurs 
ont des valeurs alEatoires, 

b) combinaison des Intensites des pixels sur lesdits 20 
supports individuellement. fournissant une intensity 
6i6mentaire, et combinaison des intensites eiemen- 
taires, fournissant une intensity filtr^e k appliquer 

au pixel courant. 

25 

Ce precede presente Tavantage de ne pas generer 
de TACHES de bruit residuel dans les images, et de ne 
pas faire apparaltre non plus des nnotifs dus h la forme 
des supports des f iltres. 

Un dispositif pour I'aide k I'observation d'une image 30 
medicale comprend : 

un systeme pour fournir des donnees d'une image 
sous forme de matrice bidimensionnelle de pixel 
ayant des valeurs d'intensite numerisees, 35 
un systeme d'affichage pour I'affichage de Umage. 
un systeme de traitement d'images ayant acces 
aux donnees dimage. et au systeme d'affichage, 
ce systeme de traitement d'image etant arrange 
pour effectuer au moins des etapes de : 4o 

a) selection, dans une fen§tre glissante cen- 
tree sur un pixel courant d un nombre (N ^ 0) 
de filtres eiementaires orientes aieatoires 
ayant des supports avec un axe d'orientation 45 
passant par le pixel courant et des caract6risti- 
ques dont une ou plusieurs ont des valeurs 
aieatoires. 

b) combinaison des intensites des pixels sur 
lesdits supports individuellement, fournissant so 
une intensite eiementaire. et combinaison des 
intensites eiementaires, fournissant une inten- 
site f iltree k appliquer au pixel courant. 

Ce dispositif presente I'avantage de restaurer une ss 
image medicale bruit6e et peu contrastee k I'origine. en 
preservant les plus fins details, et sans TACHES de 
bruit. 

L'invention est decrite ci-apres en detail en refe- 



rence avec les figures schematiques annex6es dont : 

la FIG.1 repr6sente un dispositif k rayons X pour la 
formation dimages radiologiques ou fluoroscopi- 
ques, 

la FIG.2 represente le passage d'une fen§tre 
WD(x,y) dans une image J bruit6e. 

- les FIG.3A k 3D montrent des supports de filtres 
pour mettre en oeuvre une methode de filtrage spa- 
tial par SELECTION-COMBINAISON, 

- les FIG.4A k 4D illustrent la combinaison des filtres 
des FIG.3A k 3D. 

- les FIG.5A et 5D montrent des supports de filtres 
FMH. 

les FIQ.6A k 6C illustrent la combinaison des filtres 
des FIQ.5A k 5D. 

- les FIG.7A et 7B illustrent une methode de determi- 
nation de recart type du bruit. 

1/ Dispositif k ravons X 

En reference avec la FIG.1. un syst6me de radio- 
graphie numerisee comprend une source 1 de rayons X 
; une table mobile 2 pour recevoir un patient ; un dispo- 
sitif intensificateur d'image 3 couple k un tube video 4, 
qui fournit des donnees d'images k un systeme 5 de 
traitement numerique d'image comprenant un micropro- 
cesseur. Ce dernier comprend plusieurs sorties dont 
une sortie 6 coupiee k un moniteur 7 pour la visualisa- 
tion des images radiographiques ou des images d'inten- 
site d'une zone observee du patient. 

Une image radiographique num6risee se presente 
sous forme d une matrice bidimensionnelle de pixels. 
L'image radiographique numeris6e peut contenir, par 
exemple. 512 x 512 ou 1024 x 1024 pixels codes sur 8 
bits ou 10 bits. Chaque pixel peut §tre ainsi affecte d'un 
parmi 256 ou 1024 niveaux d'intensite. Par exemple les 
regions sombres ont un niveau de gris faible, et les 
regions claires de {'image ont un niveau de gris eieve. 

L'image num6risee peut etre obtenue en mode 
fluoroscopique de rayons X. L'invention peut trouver 
une application particulierement avantageuse dans le 
filtrage d'images angiographiques. Les prises de vue 
peuvent etre effectuees sous forme de sequence d'ima- 
ges k un taux de plusieurs images par secondes, par 
exemple 30 images par seconde en mode fluoroscopi- 
que standard, ou 12.5 images par seconde en mode 
fluoroscopique puls6. 

Ueiimination du bruit sans aeation de TACHES DE 
BRUIT doit alors pouvoir Stre faite sur chaque image de 
la sequence arrivant en temps r6el. Le proc6d6 de fil- 
trage doit done 6tre rapide. et ne pas d6truire les 
details. 

Dans tous les cas. la presente invention ne tient 
pas compte de la methode par laquelle Umage numeri- 
que a ete obtenue, ni de la nature des objets qu'elle 
represente, mais concerne seulement le filtrage de 
cette image consideree seule et non en relation avec 
d'autres images de la sequence pour eiiminer le bruit 
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sans g6n6rer de TACHES de bruit 
11/ Proc6d6 de f iltraoe spatial 

En r§f6rence avec la FIG.2, le proc6d6 de filtrage 
spatial comprend la definition, dans une image J k trai- 
ter, d'un pixel appel6 pixel courant Po(xo'yo) coor- 
donn^es Xo,yo dans la matrice bidimensionnelle de 
I'image J. L1ntensit6 io(xo.yo) c*e ce pixel est suppos6e 
bruit^e au depart. 

Pour un pixel courant Po(^o.yo) d1ntensit6 Iq 
(xo.yo) dans I'image J, on d^finlt, dans une fen§tre glis- 
sante WD(xo,yo) centr6e sur ce pixel courant, une s6rle 
form6es d'un nombre N de filtres elemental res. 

Chacun des filtres de la s^rie est un filtre orients. 
Par filtre orients, on entend qu'une dimension de son 
support est plus grande que Tautre et que cette grande 
dimension, ou longueur, est parall^le k un axe appel§ 
axe d'orientation, passant par le pixel courant 
Po(xo,yo). I'autre dimension, ou largeur 6tant perpendi- 
culaire i cet axe. L'angle form6 par I'axe d'orientation 
avec, par exemple, la direction OX des lignes de la 
matrice bidimensionnelle de Timage, est appel6 angle 
d'orientation 0 du filtre el^mentaire. 

Ainsi. la serie des filtres ^l^mentaires li^s au pixel 
courant dans la fen§tre WD(xoyo) comprend plusieurs 
caracteristiques parmi lesquelles. pour chaque filtre : 

son angle d'orientation 0, 
I'etendue de son support 6gale au nombre de pixels 
paralieiemertt k I'axe d'orientation multipli§ par le 
nombre de pixels perpendiculairement k cet axe, 

et pour la s^rie de filtres : 

le nombre N des filtres, 

Chaque filtre de ta s^rie de filtres 6l6mentaires 
dans la fen§tre WD(Xo,yo) ''^e au pixel courant 
Po(xo.yo). comprend au rrioins une caract^ristTque de 
valeur al6atolre. On peut ainsl d6finlr, dans la fen§tre 
WD(xo.yo) centr6e sur le pixel courant Po(xo,yo). une 
s^rie de N filtres orient6s, dont les angles d'orientation 
0 sont al^atoires. les autres caracteristiques etant fixes. 
On peut aussi d6finir une s6rie de N filtres orient6s, 
dont k la fois les angles d'orientation © sont al6atoires 
et les etendues des supports sont al^atoires. 

En outre,^ dans la s6rie, ie nombre N de filtres peut 
etre ai^atoire. 

On d^crit ci-apr^s plus pr^cisement les caracteristi- 
ques et les proprietes du proc^de de filtrage spatial 
dans tequel on definit diff^rentes series de tels filtres 
eiementaires aieatoires. 

11,1/ Proc6d6 de sfelection-combinaison 

Dans ce proc6de app6l6 de s6lection-combinaison, 
on definit d'abord la structure des filtres ei6mentaires 
orientes ayant une ou des caracteristiques aieatoires. 
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Puis on realise une selection des pixels k Tinterieur des 
supports en fonction de leur intensite, pour seiectionner 
les pixels appartenant potentiellement au m§me objet 
que le pixel courant. Ensuite, on effectue une combinai- 

5 son lineaire ou moyenne des intensites sur chacun des 
supports des filtres eiementaires, et on combine les 
resultats de ces moyennes d'intensite, de maniere k 
fournir une intensite resultante, qui est I'intensite filtree 
k attribuer au pixel courant. 

10 Dans ce procede de seiection-combinaison, les 
supports des filtres incluent le pixel courant Pol^o.yo)- 
En reference avec la FIG.3A, on a represente, k 
titre d'exemple, la definition dune serie de N filtres 
orientes F^, F2. F3, F4, dont, pour chaque filtre : 

15 

Tangle d'orientation ©i, ©2, ©3. ©4 est fixe, 
retendue du support est aieatoire, 
les supports de ces filtres sont symetriques par rap- 
port au pixel courant PoC'^'yo)* 

20 

et le nombre N de filtres est fixe et egat k 4. 
En reference avec la FIG.36, on a represente. k 
titre d'exemple, la definition d'une serie de N filtres 
orientes F^, Fg, F3, F4,, dont, pour chaque filtre : 

25 

Tangle d'orientation By, ©2. ©3. ©4 est aieatoire. 
l etendue du support est aieatoire. 
les supports sont dissymetriques par rapport au 
pixel courant Po(xo.yo). et le nombre N de filtres 
30 est fixe et egal k 4. 

Le nombre N des filtres eiementaires peut egale- 
ment §tre aieatoire. 

En reference avec la FIG.3C, on a represente une 
35 serie de N = 3 filtres eiementaires orientes F1. F2, F3. 
dont pour chaque filtre : 

Tangle d'orientation ©1, ©2, ©3 est aieatoire, 
Tetendue des supports est aieatoire, et les supports 
40 de ces filtres sont symetriques par rapport au pixel 
courant Po(xo.yo)- 

En reference avec la FIG.3C. on peut se fixer pour 
les filtres, un nombre t^m&x d'axes d'orientation poten- 
45 tiels fixes, regulierement repartis dans Tespace autour 
de Po(xo.yo) choisissant un nombre N de fBtres 
inferieur ou egal k t^max* '^^ filtres se repartissent de 
maniere aieatoire le long des t^max d'orientation. 
Par exemple, on peut choisir N^^ax = S I ^^^^ ce cas, les 
50 8 axes d'orientation PoZ^, P0Z2. Po^a. " ^o^b po*©"- 
tiels possibles sont regulierement espacees de n/B ; et 
on peut choisir le nombre N de filtres fixe et egal k 3. 
Ces 3 filtres se repartissent alors de maniere aieatoire 
le long des 8 directions fixes possibles. 
55 Dans un autre exemple, on peut choisir le nombre 
N de filtres fixe et egal k 4. Ces filtres se repartissant de 
fagon aieatoire sur les 8 axes d'orientation. 

Dans ce type d'exemple illustre par la F1G.3C. ou 
des axes d'orientation possibles sont fixes, le nombre N 
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desf litres peut 6galement §tre al^atoire, en restant inf6- 
rieur au nombre N^ax d'axes d'orientation. 

Dans chacune des FIG.3A et 3B, I'^tendue des 
supports des filtres est de 1 x K pixels, ou K est le nom- 
bre de pixels dans le support, ce nombre K 6tant al6a- 5 
toire dans chacun de ces deux exemples. Les filtres 
repr6sent6s d titre d'exemple ont done des supports 
monodlmensionnels. Dans d'autres exemples non 
repr6sent6s. les supports des filtres peuvent §tre bidi- 
mensionnels k condition qu'une des dimensions soit au io 
moins double de I'autre dimension, cette derni^re 6tant 
petite, par exemple 2 ou 3 pixels. De pr6f6rence, la 
fenetre WD(xo,yo) est cholsie de I'ordre de 20 x 20 
pixels, ce qui fixe la grande dimension des filtres k au 
plus 20 pixels dans ce cas. is 

Parmi les caract6ristiques al6atoires potentielle- 
ment utilisabies dans la definition de la s6rie des N fil- 
tres 6l6mentaires orient§s pour atteindre aux buts de 
rinvention. c'est-a-dire pour r^aliser un filtrage spatial 
du bruit sans generation de TACHES (PATCHES) et so 
sans deterioration des details, contours et contraste, on 
propose dans l ordre decroissant d'importance : 

1) Les support des N filtres Fi d sont d'etendue 
aieatoire par exemple avec les orientations 0^ k 0^ 25 
fixes et les autres caracteristiques fixes, comme 
illustre par la FIG.3A : 

2) les supports des N filtres F^ k F^, sont d'etendue 
et d'orientation 0^ k 02 al6atoires. les autres carac- 
teristiques etant fixes, comme illustr^ par la FIG.3B 30 

3) les supports des N filtres F^ k F^ sont d'etendue 
et d'orientation 0^ k 0^ aieatoire, et en outre le 
nombre N des filtres est aieatoire, comme illustre 
par la FIG.3C ; 35 

4) les supports sont en outre asymetriques, comme 
illustr^par la FIG.3B ; 

5) un sous-6chantiIlonnage dans la fen§tre 
WD(xo,yo) peut §tre effectue prealablement au trai- 
tement de cheque pixel et le pixel de depart du 4o 
sous-6chantillonnage peut Stre aieatoire ; 

6) la moyenne MEAN pour obtenir soit les r^sultats 
individuels R^, Rg. Rn des filtres, soit le r6sultat 
general R peut s'appliquer sur des intensites pon- * 
derees par des poids al6atoires ; k cet effet, on def i- 45 
nit un coefficient aieatoire C lie au support, egal k 
rirrverse de la somme des poids, compris entre 0 et 

1 . affectes k chaque pixel dans le support ; chacun 
de ces poids peut etre egal ci 1 ; dans ce cas le 
coefficient est egal k Tinverse 1/K du nombre K de so 
pixels contenus dans le support ; le coefficient C de 
chaque filtre est pris en compte pour effectuer la 
combinaison lineaire ponderee des intensites. 

Par exemple le resultat Ri peut gtre obtenu en ss 
appliquant sur les diff6rentes intensite 1^ (xq+u, yo+v) 
du filtre FI des poids W, aieatoires et en faisant la 
moyenne MEAN de ces intensites ponder^es. d'ou il 
resulte : 



R 1 = C X L [W i(u.v) x I i(xo+u. yo+v)] 

Toutes les caracteristiques citees plus haut peuvent 
etre rendues aieatoires, seules ou en combinaison, en 
sachant neanmoins que les caracteristiques les plus 
importantes sont. pour les supports de filtres retendue 
aieatoire et I'orientation aieatoire etpour la combinaison 
le nombre de filtres N aieatoire. 

Exemple I : Le Droc6d6 de seiection-combinalson 
SIG-MEAN-MEAN 

Ce premier precede de SELECTION-COMBINAI- 
SON est d§sign6 ci-apres par le terme SIG-MEAN- 
MEAN pour sinplifier l ecriture. 

Ce procede comprend d'abord, comme d§crit d- 
dessus, la definition des filtres orientes avec : 

leurs caracteristiques aieatoires, 
leurs caracteristiques fixes, 
leur nombre fixe ou aieatoire. 

Ce procede comprend ensuite la SELECTION dans 
chaque support de filtre d'un certain nombre de pixels 
en fonction de leur difference d'intensite vis-^-vis de 
I'lntensite du pixel courant. 

Ce precede comprend ensuite la combinaison des 
filtres. 

En posant que le pixel courant Pq a pour coordon- 
n§es (xo,yo) dans la matrice bidimensionnelle de 
I'lmage J. et que le pixel P que Ton examine a pour coor- 
donnees (xq+u. yo+v), ou u et v sont des nombres 
entiers, la SELECTION est effectuee pour tout pixel P 
ayant une intensite l(xo+u, yo+v) telle que : 

IKxq+u. yo+v)-lo(^o.yo)l <mxaB (1) 

ou m est un facteur multiplicatif et ou est I'^cart type 
du bruit dans I'image (noise standard deviation, en 
anglais). 

L'ecart type du bruit dans Timage peut etre estime 
par toute methode cormue de retat de la technique, ou 
bien peut etre estim6 k priori. Uecart type du bruit est 
mesure en niveaux de gris (ou intensite). II peut §tre par 
exennple 5 ou 10 niveaux de gris. II est indispensable de 
realiser une estimation de recart type du bruit pour 
mener k bien le proc6de selon llnvention. 

On rappelle ci-aprds la definition de recart type du 
bruit telle que connue de I'homme du metier. 

L'ecart type du bruit note est l'ecart moyen que 
fait le bruit par rapport k sa moyenne. Udcart type du 
bruit est la racine carree de la variance al du bruit 
dans rimage. L'ecart type du bruit peut etre obtenu en 
appliquant la methode fournie ci-apres k titre d'exemple. 

En reference avec le diagramme de Woes de la 
FIG.7A, cette methode comprend d'abord le traitement 
de I'image J par ligne. pixel par pixel, selon un balayage 
qui peut §tre un balayage standard, depuis en haut k 
gauche de I'image jusqu'^ en bas k droite de I'image. 
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Dans une premiere ^tape symbolis^e par fe bloc MEM 
100, les intensit6s l(x,y) des pixels sont m6moris6es en 
fonction de leurs coordonn6es (x,y) dans la matrice bidi- 
mensionnelle de I'image. 

Dans une seconde 6tape symbolis^e par les blocs 
1 10Aet 1 10B, et illustr6e par la FIG7B. en chaque pixel 
courarrt P(x,y) de coordonn6es x.y d'une ligne L de 
pixels dans Timage J. on calcule deux moyennes des 
intensit^s des pixels situ^s ^ gauche du pixel courant, si 
le sens de balayage de la ligne est vers la droite : 

une premiere moyenne Ml des intensit6s sur un 
grand support Si, calcul symbolist par le bloc 
1 10A, c'est-^-dire une moyenne des Intensit6s d'un 
grand nombre de pixels de la ligne situ6s gauche 
du pixel courant. ces pixels ayant d§j^ ^t^ balay^s ; 
par grand support, oh entend environ 10 pixels, 
une seconde moyenne M2 des intensit^s sur un 
petit support S2, calcul symbolist par le bloc 11 OB, 
c'est-^-dire une moyenne des intensit6s d'un plus 
petit nombre-de pixels que celui du grand support, 
par exemple la moiti^. soit 5 pixels, et k gauche du 
pixel courant O(x.y) sur la ligne. 

Ces premiere et seconde moyennes d'intensites 
Mi. M2 peuvent §tre exactes ou pond6r§es. comme 
connu de I'homme du metier. Par pond6r6es, on entend 
qu'on affecte. k chaque intensity des pixels dont on fait 
la moyenne, un poids compris entre 0 et 1 . 

Dans une troisi^me 6tape symbolis6e par le bloc 
120 A, la m6thode de determination de T^cart type du 
bruit og, comprend ensuite le calcul d'une valeur A1 
6gale k la valeur absolue de la difference entre la pre- 
miere moyenne d'intensite M^ sur un grand support et la 
seconde moyenne d'intensite M2 sur un petit support. 
L'evaluation de cette difference : 

Ai = |Mi-M2| 

permet de rejeter les intervalles de pixels, sur' la ligne 
traitee qui contiennent des bords d'objet, c*est-^-dire 
dans lesquels apparalt une brusque variation d'intensite 
due aux contours d'un objet : on declare qu1l y a un bord 
d'objet dans I'intervalle de pixels constitue par le grand 
support S-i, lorsque la seconde moyenne d'intensite M2 
sur le petit support est tr^s differente de la premiere 
moyenne M^ d'intensite sur le grand support ; c'est e- 
dire quand ,1^ difference calcuiee entre ces deux 
nnoyennes d'intensite excede une valeur estimee en 
niveaux de gris egale k environ 10 % des niveaux de 
gris utilises pour construire I'image J numerique ; et, si 
un bord est detecte, on rejette Tintervalle constitue par 
les pixels du grand support SI . 

Dans une quatrieme etape symbolisee par le bloc 
120B. cette methode de determination de recart type 
comprend en outre le calcul d'une valeur A2 egale k la 
valeur absolue de la difference entre la premiere 
moyenne d'intensite sur le grand support et I'intensite 
du pixel courant l(x,y) 



A2 = |Mi-l(x.y)| 

L'evaluation de cette difference permet de rejeter 
les pixels aberrants. On declare que le pixel courant a 

5 une intensite aberrante lorsque la difference d'intensite 
A2 entre Tintensite du pixel courant l(x,y) et la premiere 
moyenne d'intensite M^ sur un grand support 
depasse une valeur estimee en niveaux de gris comme 
definte plus haut, environ 10 % des niveaux de gris de 

10 I'image J. 

Dans une etape ulterieure symbolisee par le bloc 
de TEST 130, cette methode comprend I'examen des 
deux conditions precedentes : 

15 si le pixel courant n'a pas une intensite aberrante, 
si I'intervalle de pixels correspondent au grand sup- 
port ne contient pas de bord d'objet. 

Dans le cas ou ces deux conditions sont conjointe- 

20 ment realisees. alors on declare que A2 qui est la valeur 
absolue de la difference d'intensite entre le pixel cou- 
rant l(x.y) et la premiere moyenne M^ d'intensite sur le 
grand support Si. est une mesure eiementaire a(x,y) 
retenue pour revaluation de recart type du bruit. 

25 La suite des operations decrites ct-dessus est repe- 
tee pour tous les pixels k traiter de Timage J. Puis une 
moyenne generale M de toutes les mesures eiemental- 
res est effectuee dans le bloc 140 et est retenue comme 
valeur de recart type du bruit os- 

30 Dans la relation (1). le facteur m multiplicatif de 
recart type du bruit os peut §tre egal e 2 ou 3. Dans le 
cas ou le bruit est gaussien et ou un filtre eiementaire 
est applique sur une zone de I'image ou I'intensite est 
sensitjiement uniforme, si m = 2, alors environ 3/4 des 

35 pixels sont seiectionnes, et si m = 3 alors pratiquement 
tous les pixels sont seiectionnes. 

On examine alors. en appliquant la relation (1). tous 
les pixels situes k Tinterieur de chaque support du filtre 
eiementaire oriente. 

40 L'operation de SELECTION d6crite ci-dessus, qui 
prend en compte recart type du bruit (lire SIGMA B) 
pour seiectionner des pixels d'intensite proche de celle 
du pixel courant, est designee ci-apres par le terme SIQ 
pour simplifier I'ecriture. 

45 En reference avec la FIG.3D, on a montre dans une 
image J, une fenetre WD(xo.yo) centree sur le pixel 
courant Po(xo.yo). des supports de f litres F^ , F2, F3. F4. 
monodimensionnels orientes dans diverses directions 
du plan de I'image, autour du pixel courant Po(xo.yo). 

50 on a mis en relief les pixels seiectionnes par SIG en les 
entourant d'un cercle. 

A Tissue de l'operation de SELECTION SIG. une 
combinaison lineaire ou moyenne des intensites des 
pixels seiectionnes est effectuee. Cette operation de 

55 moyennage est designee ci-apres par le terme MEAN. 
En reference avec la FIG.4A qui est un diagramme 
du filtrage par SELECTION-COMBINAISON SIG- 
MEAN-MEAN d'une serie de filtres telle qu'illustree par 
la FIG.3D, le filtre F^ effectue la selection des pixels 
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symbolis6es par le bloc SIG 10.1. et calcule dans le 
bloc MEAN 11.1 une premiere moyenne Ri des diver- 
ses intensit6s 

h' (xo+u, Vo+v) des pixels s6lectionn6s dans avec : 
=MEAN Il^'(Xo+u»yo+v)] 



tel que 



1 1 1 V o+u. y o+v) - 1 o(x o. y o) I < *^ X ^ B (1 ) 
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(xq+u. yo+v) d^crivant le support orient^ du liltre 

Puis, le filtre Fg calcule dans le bloc MEAN 11.2 
une seconde moyenne R2 des diverses intenslt§s 
l2'(xc+u yo+v) des pixels s6lectionn6s dans un bloc 
SIG 10 2 avec 

R „ e MEAN(l-'(xo+u, yo+v)l 

ou l intcT&iic 00c ptjictc sc*ectionn6s ob6it k la relation 
(1). (xq*^ yr.*v) deoi/ant le support orients du filtre 
F2. 

On opef e de m^e lusqu au filtre Fn compris (SIG 
10-N. MEAN 11 N) 

R^ = MEANtl^'(Xo+U,yo+v)l 

Oil les diverses mtensites In{xo+u. Vo+v) des pixels 
s6Iectionn6s obeissent a la relation (1), (xq+u, yo+v) 
d6crivant le support orlente du filtre Fn. 

Lorsque toules les moyennes Ri k Rn ont 6t6 obte- 
nues par les liltres orient§s F^ k Fn, une moyenne 
MEAN est effectuee dans MEAN 12 sur lesdites moyen- 
nes Ri ^ Rn. ce qui fournit I'intensit^ filtr6e Ro(xo,yo) 
r^ultante ^ appliquer finalement au pixel courant 

Po(xo.yo) 

Ro('^o.yo) = '^^AN [RvB2-.*Rn1 

Exemole M ■ * ^ prQC6d6 selection coifibinaison 
filG-MEAN-MED 

Ce second proc6d§ de s6lection-combinaison est 
d6sign6 ci-apr^ pour la simplicity d'§criture par le 
terme SIG-MEAN-MED. A Tissue de I'op^ration de d6fi- 
nrtion de la s6rie de N f iltres et^mentaires orient6s avec 
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intensit^s respectives lV{xo+u. yo+v) ^ In'(xo+u, yo+v) 
des f iltres Fi & Fjsj, puis les intensity Ri k Rn resultant 
de ces moyennages sont combin6es par un filtre 
median d§sign6 par MED 1 3, qui fournit I'intensit6 r6sul- 
tante Ro(xo.yo) telle que: 

Ro(^o>yo) = MED [Ri.R2 .Rn] 

Un filtre median est un filtre qui effectue le classe- 
ment des intensitds qui lui sont pr6sent6es. avec un 
num6ro d'ordre. de la plus petite d la plus grande inten- 
sity, et qui fournit comme r^sultat, rintensit§ dont la 
valeur est celle du num^ro situ6 au milieu de rensenrtble 
des num^ros, ou si les num6ros ne sont pas impairs, la 
moyenne d intensit6 correspondant aux deux num6ros 
du centre. 

Lavantage de la m6thode SIG-MEAN-MED est que 
le filtre median utilise a des propri6t6s d'injection du 
bruit haute frequence qui le rendent sp6cialemeni int^- 
ressant pour 6viter I'effet de TACHES (PATCHES) dans 
I'image. En effet le filtre median effectue la selection 
d'un des r6sultats, entre RI et RN qui est la valeur 
situ6e au milieu de ces valeurs d'intensit^, au lieu d'une 
moyenne de ces intensit^s. Le r^sultat Ro(xo.yo) Q"' 
est obtenu avec le MEDIAN est une intensity plus forte- 
ment dycorr616e, d'un pixel courant traits, au pixel cou- 
rant suivant trait^s. que rintensrt6 resultant de la 
moyenne MEAN 12 utilis6e dans Texemple I d6crit pr6- 
c§demment. 
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Exemole IM : Le proc^dy de selec tion comblnaison 
KNN 
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leurs caract6risttques al6atoires, 
leurs caract6ristiques fixes, 
. leur nombre fixe ou ai^atoire. rop6ration de selec- 
tion des pixels dans chaque support de filtre est 
effectu6e comme dans Texemple I par le crit^re (1) 
portant sur rintensit6 des pixels et prenant en 
compte r6cart type du bruit ob- 

Ensuite. en r6f6rence avec la FIG.4B, I'obtention 
des r6sultats Ri k Rn particuliers des f iltres Fi k Fn est 
effectu6e par un moyennage MEAN sur les diverses 



Ce troisi^me proc6de des selection combinaison 
est d6signy ci-apr^s, pour la simplicity d'^criture par le 
terme KNN. 

Dans cet exemple. une autre m^thode de SELEC- 
TION en fonction de intensity est appliquye au lieu du 
crityre (1) qui tient compte de rycart type du bruit oq. 
Cette autre mythode de selection consiste k retenir tous 
les pixels dont I'intensity diffyre le moins de celle du 
pixel courant. A cet effet, on calcule les diffyrences 
d'intensity 

45 A{xo.yo) " iK>^o+"- yo+v) • "o(^o'yo)i (2) 

entre les pixels dans chacun des filtres et I'intensrty du 
pixel courant. 

Dans cette mythode appeiye des K PLUS PRO- 
50 CHES VOISINS et dysignye par le terme KNN. on 
syiectionne en fait un nombre fixe K de pixels ayant 
I'intensity la plus proche de ceile du pixel courant. 

Toutes les autres caractyristiques et ytapes du pro- 
cydy KNN peuvent ytre identiques k celles de I'exemple 
55 I ou I'exemple II, et sont illustryes k cet effet par les 
Fia4C et 4D. 

Dans ces FIGS, les blocs KNN 20.1 k KNN 20.N 
symbolisent les ytapes de syiection des pixels, et les 
blocs MEAN 21 .1 MEAN 21 .N symbolisent la ryalisa- 
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tion des moyennes d'intensit^ des pixels s6lectionn6s 
sur les f litres. 

Les blocs MEAN 22.1, et MED 22.2. effectuent res- 
pectivement une moyenne ou un filtrage median des 
r^sultats particuliers des fillres §l6mentaires pour four- 5 
nir le r^sultat final Ro(xoyo)- 

11.2/ Le proc6d^ FMH 

Dans les 6tapes du proc6d6 d^crit ci-apr6s on d6fi- io 
nit une s6rie de N filtres Iin6aires 6l§mentaires moyen- 
neurs. A r^ponse impulsionnelle, dont les r^ponses sont 
en outre combin6es par des filtres medians, ce qui 
revient k designer cette classe de filtres par FMH (en 
anglais : FIR HYBRID MEDIAN FILTERS). Cette combi- is 
naison a une structure d'arbre qui comprend les filtres 
lin^aires moyenneurs (FIR) et les filtres medians 
(MED). 

En r6f6rence avec les FIG.5A d 5D. on d^finit une 
fen§tre WD(xo,yo) bidimensionnelle, centr^e sur le so 
pixel courant Po(xo-yo) coordonn6es {xQ.yo) dans la 
matrice bidimensionnelle de Timage J k trarter. Dans 
cette fen§tre WD(xo.yo). on d^f init des supports de fil- 
tres selon les lignes des pixels, et les colonnes de 
pixels, et les deux diagonales entre les lignes et les 25 
colonnes. Le pixel central Po(xo.yo) de la fen§tre 
WD(xo.yo). Q»J' est le pixel courant. est consld6r§ ^ 
part. 

-En r6f §rence avec la FIG.5A, on a repr§sent6 k titre 
d'exemple une s6rle de N filtres orient6s H^ H2, V^, Vg, 30 
' ^11. ^^2' D21. O22 dont pour chaque f litre : 

I'orientatlon est fixe consid4r6e par exemple par 
rapport k Taxe OX, selon la verticale e^, I'horizon- 
tale 02 (angle non repr6sent6) et les diagonales 35 
©3. ®4 passant par le pixel courant. centre du carr6 
formant la fen§tre, 

r^tendue du support est al6atoire : ces supports 
sont monodimensionnels. c*est-&-dire ont une lar- 
geur de 1 pixel et une longueur qui est Ilmit6e par le 40 
cadre de la fenStre WD (xo,yo). 
les supports des filtres H^. H2 disposes selon la 
m§me direction, par exemple horizontale, sent 
sym^triques. II en est de m§me des supports de fil- • 
tres V-i et V2 disposes selon la direction verticale, 45 
des supports de filtres D^^ et D^g disposes selon la 
premiere diagonale, des supports de filtres D12, 
D22 disposes selon la seconde diagonale, 
le nombre N des filtres est fixe ; Ici ce nombre est 8, 
et ces supports de filtres sont disposes selon tous 50 
les angles possibles d'orientation 0^ k 04 multiples 
de rc/4 dans I'espace. 

dans la s6rie de filtres. on incorpore le pixel central 
Po(xo.yo) de la fenStre et son Intensity est consid^- 
r6e k part. 55 

En r6f6rence avec laFIG.SB, une s6rle de filtres du 
mime type est representee avec les caract^rlstiques 
suivantes, par rapport k la FIG.5A : 



• les supports des filtres qui sont disposes selon la 
mime direction sont dissym^trlques de mani^re 
aieatolre. 

En reference avec les FIG.5C el 50, au lieu de con- 
siderer un nombre de filtres disposes selon tous les 
angles potentiels d'orientation 0^, 03, ©3, 04, seule- 
ment deux axes sont considerees : les axes horizontal 
et vertical, sur la FIG.5C ; les axes formes par les direc- 
tions des deux diagonales sur ta FIG.5D. D'ou, sur la 
FIG.5C, la s6rie de filtres comprerxl seulement les sup- 
ports Hi, H2, V^, V2 alignes sur ces axes horizontal et 
vertical ; et sur la FIG.5D. la serie de filtres comprend 
seulement les supports D^^, D^g- ^21 ^22 a'ignes sur 
ces axes en diagonale. Sur ces FIG.5C et 5D k titre 
d'exemple, les supports ont une etendue al6atoire et 
sont symetriques par rapport au pixel courant. 

En reference avec les FIG.6A k 6C, les filtres ayant 
les supports de caract6ristiques al6atoires decrits en 
relation avec les FIG.5A k SO effectuent d'abord la com- 
binaison Iin6aire ou moyenne MEAN des intensltes des 
pixels sur teur support. 

On a appeie l'vi. I'v2. I'hi. I'h2 '©s differentes inten- 
sltes des pixels situes sur les supports V^, Vg, H^ et H2 
respectivement. 

On a appeie I'dh. I*di2. I'dsl '*d22 les differentes 
intensites des pixels sltu6s sur les supports D11. D12, 
D21, D22 respectivement. 

Les resultats des moyennes effectuees, dans les 
blocs MEAN31.1 k MEAN 31.8. par les filtres H^. H2, 
Vi. V2, D11. D12. D21. D22 sont des intensltes notees 
'hv 'h2' 'vi. 'v2. 'dii' 'di2. >D21 ®* Id22 respectivement. 

En reference avec la FIG.6A. la combinaison en 
arbre des resultats des filtres nnoyenneurs MEAN pre- 
cedents comprend 2 niveaux incluant I'introductlon 
dans un filtre median MED 32.5 k 5 valeurs, des inten- 
sltes moyennes Ivi, Iv2. 'hi» Ih2 sur les axes horizontal 
et vertical plus I'intenslte de Iq (xq. yo) du pixel central ; 
et en paralieie, I'introduction dans un filtre median MED 
32.6 k 5 valeurs. des intensites nroyennes Iqh. bij. 
Id2i ' '022 sur les deux directions diagonales plus rinten- 
site lo{xo, yo) du pixel central. Ces deux filtres medians 
fournissent respectivement les intensites resultantes 
R1 el R2, termlnant le premier niveau de I'arbre. 

Les deux Intensltes R1 el R2 calcul6es au premier 
niveau de Tarbre sont introduites. avec I'intenslte {©(xq. 
yo) du pixel central dans un nouveau filtre median MED 
33, k 3 valeurs, qui fournit Tintensite resultante 
Ro{xo.yo) ^ appliquer comme valeur d'intensite filtree 
sur le pixel courant Pq (xq. y©)- 

En reference avec la FIG.6B. la combinaison com- 
prend egalement 2 niveaux, mals elle n'est appllquee 
que sur les moyennes d'intensite Ivi, Iv2. 'hi. *H2 relati- 
ves aux seules directions vertlcales et horlzontales res- 
pectivement. Cette combinaison comprend 
introduction des moyennes d'intensite Ivi. Iva plus 
I'intenslte lo(xo.yo) brultee du pixel central dans un filtre 
median MED 32.1 a 3 valeurs pour fournir une Intensite 
resultante R1 ; et I'introductlon des moyennes d'lnten- 
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sit6 Ihi . Ih2 en parallfele dans un filtre median MED 32.2 
a 3 valeurs pour fournir une intensity r^sultante R2. 

Les intensit6s R1 el R2 calcul^es au premier 
niveau de I'arbre sont introduites, avec I'lntensit^ lo(xo. 
Yo) du pixel central dans un nouveau filtre median MED s 
33.1. k 3 valeurs. qui fournit I'intensit^ r^sultante 
Ro(xo.yo) ^ appliquer comme valeur d'intensit6 filtrde 
sur le pixel courant Pq (xq. Vo)- ^• 

On peut. dans une combinaison non representee, 
renplacer les valeurs de moyennes d'intensite Ivi. Iv2. 
Ihi. Ih2 dans la FIG.6B par les moyennes obtenues 
dans les m§mes conditions sur les deux diagonales, 
soit respectivement Iqh, Idi2' b2i» 'd22' POur former 
une autre combinaison de flltres & deux niveaux sur 
seulement deux axes qui sont celies des deux diagona- is 
les. 

En reference avec la FIG.6C. la combinaison en 
arbre comprend 3 niveaux. Les deux premiers niveaux 
sont constitu6s comme d^crit en relation avec la FIG.6B 
de la combinaison des filtres sur les deux axes horizon- so 
tal et vertical et donne par MED 33.1 un r^sultat R3 ; et 
en paralieie comme d6crit plus haut. de la combinaison 
des filtres sur les deux axes diagonaux et donne par 
MED 33-2 un r^sultat R4. 

Le premier r^sultat R3 et le second r^sultat R4 sont 2S 3. 
ensuite combines avec I'lntensite bruit^e Iq (xq. yo) 
pixel courant dans un filtre median MED 34 ^ 3 valeurs. 
pour donner en sortie de ce troisi^me niveau, I'intensite 
r6sultante Ro(xo.yo) ^ appliquer comme valeur d'inten- 
sit6 filtr^e sur le pixel courant Pq (xq. yo)- 

Lintensite d un pixel courant est f iltr6e par la s6rie 4. 
de filtres ei^mentaires al6atoires pr6alablement d6finis. 
selon les 6tapes de selection et de combinaison propo- 
s6es. puis les intensit6s des autres pixels d6sir6s de 
I'rmage sont f iltr^es, jusqu'^ ce que tous les pixels d6si- 35 
r^s de I'image aient subi le traitement selon le proc^de 
propose, c'est'^-dire. en operant le f iltrage de Tintensite 
de chaque pixel desire par une serie de filtres 6l6men- 5. 
taires dont les caracteristiques qui sont fixes ou qui sont 
aieatoires ont 6te prealablement definies, ainsi que le 40 
mode de realisation des combinaisons d'intensite. 

Revendt cations 

1 . Precede de f iltrage spatial du bruit dans une image 45 
(J) ayant la forme d'une matrice bidimensionnelle 
de pixels [P(x,y)] ayant des valeurs d'intensite 
[|(x,y)] numerisees bruit6es. ce procede compre- 
nant des etapes de : 

so 

a) selection, dans une fenetre glissante 
[WD(xo. yo)] centree sur un pixel courant 6. 
■ [Po(xo. yo)] d'un nombre (N ^ 0) de filtres 616- 
mentaires [F^, F2.... F^ : V^. Vg. Hi, H2. D^i, 
Di2. D21, D22] orientes aieatoires ayant des ss 
supports avec un axe d'orientation [PqZi, 

P0Z2 Pq^n] passant par le pixel courant et 

des caracteristiques dont une ou plusieurs ont 
des valeurs aieatoires. 



b) combinaison des intensites des pixels sur 
lesdits supports individuellement. fournissant 
une intensite eiementaire. et combinaison des 
intensites eiementaires. fournissant une inten- 
site f iltree [Ro(Xo. yo)] ^ appliquer au pixel cou- 
rant, 

Procede selon la revendication 1 , selon lequel les 
supports des filtres ont potentiellement. parmi les 
caracteristiques dont une ou plusieurs ont des 
valeurs aieatoires : 

teur nombre (N) dans la serie, 

un angle d'orientation [©i. 02-" ®n1 entreTaxe 

d'orientation [PqZi. P0Z2 Pq^nI ^^e 

de reference (OX) defini dans la matrice bidi- 
rectionnelle, 

une etendue definie en nombre de pixels dont 
la plus grande dimension est paralieie k Taxe 
d'orientation, 

une asymetrie par rapport au pixel courant, 
un coefficient lie k des coefficients de ponddra- 
tion appliques aux intensites des pixels. 

Precede selon Tune des revendications 1 ou 2, 
selon lequel un sous-echantillonnage des pixels de 
la fenetre est potentiellement pratique avec un pixel 
de depart de ce sous-echantillonnage localise de 
maniere aieatoire. 

Precede selon Tune des revendications 1 k 3, selon 
lequel lesdits supports de filtre sont monodimen- 
sionnels avec les pixels alignes sur les axes 
d'orientation, ou sont bidimensionnels avec la 
dimension du support paralieie k I'axe d'orientation 
au moins double de la dimension transverse. 

Precede selon Tune des revendications 1 k 4, selon 
lequel : 

les angles d'orientation potentiels [©^ ©2. .. 
0m] sont fixes et regulierement r6partis dans 
I'espace, 

le nombre (N), fixe ou aieatoire. des filtres eie- 
mentaires aieatoires est inferieur ou egal au 
nombre (Njnax).des angles d'orientation, 
et les axes d'orientation des supports de ces fil- 
tres eiementaires sont disposes selon ces 
angles d'orientation de maniere fixe ou aiea- 
toire. 

Precede selon la revendication 5, selon lequel : 

les supports des filtres eiementaires aieatoi- 
res incluent le pixel courant. 
ce precede comprenant en outre, entre retape a) et 
retape b), une etape de : 

a') selection (SIG. KNN) des pixels en fonction 
de leur intensite sur chaque support de filtre 
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6l6mentaire. selon un crit^re qui retient les 
pixels d'intensit§ proche de celle du pixel cou- 
rant pour estimer leur appartenance k un 
m§me objet de Timage. 



Iimage. 

14, Disposltif pour I'aide d robservatton d'une image 
m^dicale, comprenarrt : 



7. Proc6d§ selon la revendication 6, selon lequel, pour 
la selection des pixels, le crit^re retient les pixels 
qui ont une intensity [l(xo+u. Yo+v)] telle que la 
valeur absolue de la difference entre ladlte intensity 

et I'intensite [lo(xo. Vo)] du pixel courant est inf6- io 
rieure ^ un multiple (m) de I'^cart type du bruit (ag) 
dans rimage (J). 

8. Procede selon la revendication 6, selon lequel : 

pour la selection des pixels, le crit^re retient is 
les pixels en un nomtKe predetermine (K) dans le 
support 

9. Procedo sckxi lunc* des revendications 7 ou 8, 
selon lepuci ta combmarson des intensit6s des so 
pixels sur les &tx)port& est hneaire et la combinai- 
son des int ensues ©lementaires est effectu6e par 

un filtre moyenneur (MEAN) 

10. Procede selon I une des revendications 7 ou 8, 25 
selon lequel la comtxnaison des intensites des 
pixels sur les supports est lin^aire et la combinai- 
son des intensrtes elementaires est effectufee par 

un filtre median (MED). 

30 

11. Procede selon la revendication 5, selon lequel : 

les axes d'orientation fixes sont au nombre de 
quatre. disposes selon la ligne, la colonne, et 
les deux diagonales passant par le pixel cou- 35 
rant, dans la fenetre, 

les supports sont positionnes pour exclure le 
pixel courant, etant disposes de part et d'autre 
du pixel courant, sur la ligne (H^, H2). sur la 
colonne (V^. V2), sur la premiere diagonale 40 

D12) et sur la seconde diagonale (D21. 
D22) passant par le pixel courant, 
le pixel courant, central dans la fenetre [Wq 
(xq, Vo)] . est considere comme un support par- • 
ticulier, 45 



un systeme pour fournir des donn6es d'une 
image sous forme de matrice bidimensionnelle 
de pixels ayant des valeurs d'intensite numeri- 
sees. 

un systeme d'affichage pour I'affichage de 
rimage, 

un systeme de traltement d'images ayant 
acces aux donnees d'Image, et au systeme 
d'affichage, 

ce systeme de traitement dimage etant arrange 
pour effectuer des etapes selon le procede con- 
forme e Tune des revendications 1 e 13. 



et la combinaison des intensites elementai- 
res prend en compte I'intensite du pixel courant. 

12. Procede selon la revendication 11, selon lequel la 
combinaison des intensites des pixels sur les sup- 
ports est lineaire. et la combinaison des Intensites 
eiementaires est effectuee par des filtres medians 
sur deux ou trois niveaux, prenant. en compte 
I'intensite du pixel courant. 

13. Procede selon Tune' des revendications 1 e 12, 
selon lequel les etapes def inies pour un pixel e trai- 
ler sont repetees pour tous les pixels k trailer dans 
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